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Abstract. The system size and centrality dependence of multiplicity fluctuations in nuclear colli- 
sions at 158A GeV as well as the energy dependence for the most central Pb + Pb collisions were 
studied by the NA49 experiment at CERN SPS. A strong increase of fluctuations was observed with 
decreasing centrality in C + C, Si + Si and Pb + Pb collisions. The string hadronic models (UrQMD, 
Venus, HIJING, HSD) can not reproduce the observed increase. This may indicate a strong mixing 
of target and projectile contribution in a broad rapidity range. For the most central collisions at all 
SPS energies multiplicity distributions are significantly narrower than a corresponding Poisson one 
both for negatively and positively charged hadrons. The UrQMD model seems to reproduce the 
measured values on scaled variance. Statistical model calculations overpredict results when conser- 
vation laws are not taken into account. 
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INTRODUCTION 

At high energy densities (~ I GeV / fm 3 ) a phase transition from hadron gas to quark- 
gluon-plasma (QGP) is expected to occur. There are indications that at top SPS energies 
quark-gluon-plasma is created at the early stage of heavy ion collisions [ 1 ] . Lattice QCD 
calculations suggest furthermore the existence of a critical point in the phase diagram 
of strongly interacting matter which separates the line of the first order phase transition 
from a crossover. Models predict an increase of multiplicity fluctuations near the onset 
of deconfmement [2] or the critical point |0]. Statistical model calculations [4] showed a 
non-trivial decrease of fluctuations due to conservation laws. These reduced fluctuations 
would then serve as a "background" in the search for the critical point. 



THE NA49 EXPERIMENT 

The NA49 detector is a large acceptance fixed target hadron spectrometer described 
in |@]. The centrality of a collision is determined using a downstream Veto calorimeter 
which measures the energy in the projectile spectator domain [6]. This allows a deter- 
mination of the number of projectile participant nucleons Np WJ . The number of target 
participants is not fixed in the experiment, model calculations show large fluctuations of 
their number in peripheral collisions 02D. 



Analysis Procedure 



Since detector effects like track reconstruction efficiency might have a large influence 
on multiplicity fluctuations, it is important to select a very clean track sample for the 
analysis. Therefore the acceptance for this analysis is limited to a part of the forward 
hemisphere, where the NA49 detector has the highest tracking efficiency. This was done 
by restricting the analysis to the rapidity interval 1 < y(n) < ybeam 1 f° r 20A to 80A 
GeV and 1.08 < y(n) < 2.57 for 158A GeV. In addition a cut on transverse momentum 
according to [SJ] was applied, which is dependent both on rapidity and azimuthal angle. 
For such tracks the reconstruction efficiency is larger than 98%. 

The basic measure of multiplicity fluctuations used in this analysis is the scaled 
variance: CO = > where Var(n) and < n > are variance and mean of the multiplicity 
distribution, respectively. In this paper only the results for negatively charged hadrons 
are shown. 

The data is corrected for the finite width of centrality bins. Data on centrality and sys- 
tem size dependence at 158A GeV is corrected for the resolution of the Veto calorimeter, 
for most central Pb + Pb data this correction is small (it would decrease the results by 
less than 5%) and was neglected. 

The total systematic error is estimated to be +2% respectively —5% [6] for the most 
central Pb + Pb collisions, the systematic error for p + p, C + C, Si + Si and non-central 
Pb + Pb collisions is shown in figured 

CENTRALITY AND SYSTEM SIZE DEPENDENCE 

The centrality dependence of the scaled variance for negatively charged hadrons (h~) 
at 158A GeV is shown in figure [U Scaled variance increases with decreasing centrality 
of a collision. For very peripheral collisions there is a hint that it decreases again, but 
the systematic errors are large in this region. Data is compared to predictions of various 
string hadronic models for Pb + Pb collisions; they all predict, in contradiction to data, 
a flat centrality dependence. Scaled variance behaves similar in p + p, C + C, Si + Si 
and Pb + Pb collisions if plotted against centrality defined as Np roj /A, where A is the 
number of nucleons in nuclei. This is a hint that not the size of the collision system 
is the correct scaling parameter for the observed increase of CO but the fraction of the 
colliding nucleons in the projectile nuclei. In [9] it is suggested that target participants 
contribute to particle production in the projectile hemisphere and their fluctuations 
therefore result in an increase of multiplicity fluctuations. The predictions of this effect 
called mixing is in approximate agreement with data as shown in figure [T] (right). Inter- 
particle correlations II 1 Oil and percolation models [1JJ provide alternative explanations 
of the results. 



Rapidity is calculated in the center of mass system assuming pion mass. 
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FIGURE 1. Left: Comparison of scaled variance of h~ in Pb + Pb collisions (circles) with string 
hadronic models: HIJING 0, UrQMD, HSD H and Venus Q. Middle: Centrality dependence of 
scaled variance of hr in p + p, C + C, Si + Si and Pb + Pb collisions at 158A GeV. Right: Scaled variance 
of h~ in comparison to model calculations for transparency, mixing and reflection of matter in the early 
stage |9]. The outer errors correspond to a sum of statistical and systematical uncertainties. 
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FIGURE 2. Energy dependence of multiplicity fluctuations of h~ in Pb + Pb collisions in comparison 
to string-hadronic models Venus 1 1 311 and UrQMD 1 14] (left) and in comparison to canonical and grand- 
canonical statistical hadron-resonance gas models 11511 (right). Only statistical errors are shown, the 
systematical errors (not shown) are aVy S %. 



ENERGY DEPENDENCE IN CENTRAL COLLISIONS 

At all energies the scaled variance of hr is significantly smaller than one, the value for 
a corresponding Poisson distribution (see figureEl) 1 6]. A direct quantitative comparison 
of CO at different energies is not possible due to different experimental acceptance. The 
UrQMD model is in approximate agreement with data. The Venus model overpredicts 
0) at higher energies. No significant increase of multiplicity fluctuations due to critical 
point or onset of deconfinement is observed. 

Predictions of a statistical model II 1 511 for canonical and grand-canonical ensembles 
are compared to data. The measured scaled variance is much lower than predicted by the 
grand-canonical ensemble. The canonical model predicts co smaller than 1, in qualitative 
agreement with data. Energy momentum conservation and the finite volume of hadrons 
are expected to cause an additional suppression of fluctuations. Thus the observed small 
fluctuations seem to be a non-trivial effect of conservation laws in a relativistic hadron 



gas i]. 
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